Linear and Nonlinear Optical Constants of BiFeOs 
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Using spectroscopic ellipsometry, the refractive index and absorption versus wavelength of the 
ferroelectric antiferromagnet Bismuth Ferrite, BiFeOa is reported. The material has a direct band- 
gap at 442 nm wavelength (2.81 eV). Using optical second harmonic generation, the nonlinear optical 
coefficients were determined to be di5/d22 = 0.20±0.01, d3i/d22 = 0.35±0.02, £^33/^22 = -11.4±0.20 
and |d22| = 298.4±6.1 pm/V at a fundamental wavelength of 800 nm. 
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BiFeOs (BFO) is an antifcrromagnetic, ferroelectric 
with Neel temperature T^v = 643 K, and ferroelectric 
Curie temperature Tc — 1103 K.P, 0, Q It is presently 
one of the most studied multiferroic materials due to its 
large ferroelectric polarization of ~ 100 ^C/crn^ in thin 
films, and the possibility of coupling between magnetic 
and ferroelectric order parameters, thus enabling manip- 
ulation of one through the other. [l| Linear and nonlinear 
optical spectroscopy tools are ideally suited to study such 
coupling. ^ While the mean refractive index for bulk sin- 
gle crystal BiFeOs has been previously investigated, 
the optical constants of thin films have not been pre- 
sented thus far. Also, an indirect gap at 673 nm (1.84 
eV) was reported beforeQ, which is shown here to be 
an absorption onset potentially due to a joint density 
of states effect and not associated with phonon partic- 
ipation. In our analysis, the material appears to have 
a direct gap with a bandedge at 442 nm instead. No 
studies of nonlinear optical coefficients of BiFeOs, in any 
form, exist. In this letter, we measure large second order 
optical nonlinearities in BFO. 

Epitaxial and phase-pure BFO thin films were syn- 
thesized by pulsed-laser deposition (PLD) as well as 
molecular-beam epitaxy (MBE) on (111) SrTiOs (STO) 
substrates. [H 0] The films studied here are epitaxial 
with orientation relationship BFO(0001)//STO(lll)and 
[2llO]BFO//[110] STO. We note specifically that unlike 
many epitaxial thin films, these films do not have any ad- 
ditional structural variants, including any rotational vari- 
ants within the film growth plane. Thus, these (0001) ori- 
ented films have nearly single crystalline perfection, with 
three well-defined crystallographic a;-[2Tl0], and ^/-[liOO] 
axes within the film plane, and the z-[0001] axis normal 
to the plane. The three y-z mirror planes in the 3to point 
group symmetry for BFO are thus well defined and al- 
low us to extract nonlinear coefficients precisely without 
ambiguity. Typical film stoichiometry, as determined by 
Rutherford backscattering spectrometry (RBS), was sto- 
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71-3 : E„ = 4.1 1 ± 0.02 eV 
r = 1 .71 ± 0.06 eV , A = 61 .74 ± 5.00 eV 



TL2 : E„ = 3.01 ± 0.01 eV 
■- r =0.57± 0.03 eV , A= 31.36 ± 2.22 eV 



TL1 : = 2.46 ± 0.01 eV 
' r =0.36 + 0.05eV, A= 11.60±2.4OeV 
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FIG. 1: (a)Index of refraction (n) and extinction coefficient 
(k) for BiFe03 deposited on (111) SrTiOs over a spectral 
range from 190 to 1670 nm. Also shown are the dielectric 
function parameters for the three Tauc-Lorentz oscillators 
(TL 1, TL 2, TL 3) within the spectral range. The param- 
eterization also includes an additional oscillator outside the 
spectral range with Eq = 7.29 ± 0.06 eV, F = 4.28 ± 0.34 eV 
and A = 51.19 ± 2.78 eV. (b) Plot of (aE)!''^ and {aEf vs. 
photon energy E where a linear extrapolation of (aE)^ to 
suggests a direct bandgap at 2.81 eV. 



ichiometric within 3 % error of the measurement (Bi:Fe 
— 0.98-0.99:1). There were no amorphous or secondary 
phases as confirmed by transmission electron microscopy. 
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(a) Tilt e = 0° P (b) Tilt B = + 45" 



FIG. 2: Variation of p and s polarized SHG intensity with 
incident polarization angle for a BF0//ST0(111) film with 
E-axis perpendicular to the plane of incidence in (a) Normal 
incidence 9 = 0° and (b) 9 = 45° tilt. 



Ellipsometric spectra in (A, were collected ex situ 
for a BiFeOa film prepared by MBE on (111) SrTiOs at 
Oi = 55° and 70° angles of incidence using a variable- 
angle rotating-compensator multichannel spectroscopic 
ellipsometerQ with a spectral range from 190 to 1670 
nm. The optical properties (n, k) shown in Fig. [TJa) and 
the corresponding dielectric function spectra (£1,62) are 
extracted by using a least squares regression analysis and 
a weighted root mean square error [9], to fit the ellipso- 
metric spectra to a four-medium optical model consist- 
ing of a semi-infinite STO substrate / bulk film / surface 
roughness / air ambient structure. The free parameters 
correspond to the bulk and surface roughness thicknesses 
of the film and a parameterization of the BiFeOa dielec- 
tric function. The dielectric function parameterization 
of BiFcOa consists of four Tauc-Lorentz oscillators [l3| 
sharing a common band gap and a constant additive 
term to ei denoted by £00 (equal to 1 for this model). 
The parameters corresponding to each oscillator include 
an oscillator amplitude A, broadening parameter F, res- 
onance energy Eq, and a Tauc gap Eg common to all os- 
cillators. The optical properties of the surface roughness 
layer are represented by a Bruggeman effective medium 
approximation consisting of a 0.50 bulk film / 0.50 void 
mixture. ^ This model yields the common Tauc gap Eg 
= 2.15±0.06 eV, bulk thickness d6= 468.93±0.78 A, and 
surface roughness thickness ds= 75.39±0.4 A. 

We note that though the Tauc gap at 2.15 eV (577 
nm) represents the onset of absorption, it is not the in- 
direct gap, as claimed in literature. @ A plot of a^E^ 
vs. photon energy E (a = 47rk/A) and the linear ex- 
trapolation to E^ = indicates a direct gap at 2.81 
eV (442 nm) as shown in Fig. [TJb). This value is in 
good agreement with that obtained from more recent 
optical measurements. 12] Band gap measurements on 
different MBE-grown BiFeOa films grown on (001) and 
(111) SrTiOa substrates revealed a direct band gap in 



all cases with Eg = 2.77±0.04 eV. The presence of two 
distinct slopes in (aE)-'^/^ vs. E characteristic of an indi- 
rect band gap [l^ is not observed. We obtain the linear 
complex indices from this model to be iV^ =2.836+0i and 
iV|"=3. 444-1-0. 981i for corresponding wavelengths of 800 
and 400 nm, respectively. It should be noted that al- 
though BiFeOa is uniaxially anisotropic, only the optical 
properties of the ordinary index of refraction have been 
obtained for this film.[l3| 

The crystal symmetry of epitaxial BFO(lll) films has 
been shown to be point group 3m using optical SHG 
and diffraction techniques. Optical SHG[i^ involves 
the conversion of light (electric field iJ") at a frequency 
Lu into an optical signal at a frequency 2lu by a nonlin- 
ear medium, through the creation of a nonlinear polar- 
ization P^'^ oc dijkEj' E'^ , where dijk represent the non- 
linear optical coefficients. BFO film with thickness of 
about 50 nm grown on STO(lll) substrates was used for 
this study. STO is centrosymmetric (cubic) and does not 
contribute SHG signals of its own for the incident powers 
used. The SHG experiment was performed with a funda- 
mental wave generated from a tunable Ti-sapphire laser 
with 65 fs pulses of wavelength 800 nm incident from 
the substrate side at variable tilt angles 9 to the sample 
surface normal. 

As shown in Fig. [21 the crystallographic y-z plane in 
the BFO film was aligned with the incidence plane. The 
polarization direction of incident light is at an angle (fi 
from the x axis, which was rotated continuously using a 
half-wave plate . The intensity J^" of the output SHG 
signal at 400 nm wavelength from the film was detected 
along either j = p, s polarization directions as a function 
of polarization angle (j) of incident light. The resulting 
polar plots of SHG intensity for p and s-polarized output 
a,t e = 0° and 9 = 45° are shown in Fig. ^a) and (b) 
respectively. If the incident beam has intensity Iq then 
the nonlinear polarizations for BFO(lll) film with x-axis 
perpendicular to the plane of incidence is given bv p^lT^ . 

Pi^^ = /o/. sin 2(/)(di5/. sin - ^22/,, cos 9) 

Py"- = h{-d22 cos^ cffl + d22fl cos2 6 sin^ 

+ (ii5/y/zsin26'sin^(/)) 

Pf ^ = /o (-^31 cos2 4>fl + dai /' cos^ 6 sin^ 

+ daa/'sin^^sin^c/,) (1) 

where d.y are nonlinear coefficients and fi are effective 
linear Fresnel coefficients. The measured intensity of the 
p and s polarized SHG in transmission geometry (ne- 
glecting birefringence) is proportional to nonlinear po- 
larization. The expected SHG intensity expressions for p 
and s output polarizations in the predicted 3m symmetry 
system of BFO are: 

ll'^ = A{cos^ (t) + B sin^ 0)2 

Csin2 20 (2) 

where B and C are given by 
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B = 
C = 



Ki5fyfz sin 29 cos 9b + K^if^ cos^ 9 sin Ob + K^^fl sin^ 6 + cos^ cos 

K:iifl sin 6*5 - cos 9b 
D{Ki5fJ^ sin 6* - f^fy cos 6*) 



(3) 
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FIG. 3: Variation of extracted B parameter from p polar- 
ized SHG signal {Bexpt = —[-fi^j^^pr^]^) and C Parameter 

{Cexpt = Il'^{<t> = 45°)) as a function of tilt angle, 0. The 
solid lines show theoretical fits to the data. 



the film, respectively, p^^^p"^) are the second- harmonic 
(fundamental) signal powers measured, T is the trans- 
mission coefficient of the fundamental, A is the area of 
the probed beam, n are the indices of refraction, 1^ the 
coherence lengths and a = 47rk/A is the absorption coef- 
ficient at 2u}. For this experiment, Af = Ar = 7r(60/im)^ 
and the coherence length of the film Zc(filni) is equal to 
the film thickness. 

The B and C parameters, which contain the ratios of 
nonlinear coefficients, are experimentally obtained by col- 
lecting the p-in-p-out /p"(0 = 90°), s-in-p-out /p'^((/' = 
0°) and 45-in-s-out = 45°) SHG signals for dif- 

ferent angles of tilt 9 about the x axis. The experimen- 
tal data for B and C parameters (Fig. [3]) is then fit- 
ted to Eq. [3] to extract the ratios K15 = 0.20 ± 0.01, 
K31 = 0.35 ± 0.02 and K33 = -11.4 ± 0.20. Taking ab- 
sorption into account, the estimated effective coefficients 
are 



Here K15 = ^15/^22, K31 = ^31/^22 and K33 = 
^33/^22 are the ratios of the nonlinear optical coefficients, 
A and D are scaling parameters, and is the angle that 
the generated second harmonic wave makes with the sur- 
face normal inside the film. Theoretical fits to the ex- 
perimental polar plots based on Eq. [2] are excellent both 
in normal incidence and tilted configuration as shown in 
Fig. [^Ka) and (b), respectively, for both p and s polar- 
ized SHG output. In normal incidence {9 = 0°), only the 
1^22 coefficient is involved in the generated Ip'^{4> = 0°). 
Using the 1^22 =1.672 pm/V coefficient of a single crystal 
z-cut LiTa03 used as a reference, the d22 coefficient of 
the BEG film is calculated by employing the following 
equation [l^, 





(4) 

where the subscripts r and / refer to the reference and 



1^22!= 298.4 ±6.1 pm/V, 1^31 1 = 59.7 ± 4.2 pm/V 
Idisl = 104.4 ± 8.1 pm/V, 1^33! = 3401 ± 129 pm/V. 

Note that only the signs of the ratios if 15, -K'31 and ^^33 
were determined unambiguously. The absolute signs of 
the dij coefficients were not determined, except to state 
that the 1^33 coefficient has the opposite sign to the other 
coefficients. The large values of dij coefficients most 
likely arise due to electronic resonances at the 400 nm 
SHG wavelength. 

To conclude, we report the complex index of refraction 
versus wavelength and optical second harmonic genera- 
tion coefficients in BiFe03 thin films. These studies will 
be important in performing further linear and nonlinear 
optical spectroscopy of the magnetism and ferroelectric- 
ity in this material. 
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